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Abstract. Recent calculations of evolutionary tracks of TP-AGB stars of different mass and 
metallicity by Marigo et al. (2007) have been incorporated in the Bruzual & Chariot evolution- 



ary population synthesis models. The mass of the stellar population in HUDF galaxies at z from 
1 to 3 determined from fits to the spectro-photometric data of these galaxies is 50 to 80% lower 
than the mass determined from the BC03 models. The ages inferred for these populations are, 
with exceptions, 40 to 60% of the BC03 estimates. 



1. Introduction 

In the oral presentation of this paper I gave a summary of the current status of various 
high resolution stellar spectral libraries and their use in population synthesis models. A 
written version of this summary can be found in a series of conference papers appeared 
in recent years, (Bruzual 2004, 2005, 2006) and will not be repeated here for reasons of 
space. In this paper I will concentrate in the second part of my talk which versed on 
the role of the thermally pulsing asymptotic giant branch (TP-AGB) phase of stellar 
evolution on the integrated properties of stellar populations. 

It has been pointed out by several authors, e.g. Maraston et al. (2006) | Kannappan & Gawiser (2007) 
that the estimates of the age and mass of the stellar population present in a galaxy de- 
pend critically on the ingredients of the stellar population model used to fit the galaxy 



spectrum. Maraston et al. (2006) have shown that the treatment of the TP-AGB phase 



of stellar evolution is a source of major discrepancy in the determination of the spec- 
troscopic age and mass of high-z (1.4 < z < 2.7) galaxies. The mid-UV spectra of these 
galaxies indicate ages in the range from 0.2-2 Gyr, at which the contribution of TP- 
AGB stars in the rest-frame near-IR sampled by Spitzer is expected to be at maximum. 



Maraston et al. (2006)| find that in general the Maraston (2005) models (M05 hereafter) 
provide better fits to the observations than the Bruzual & Chariot (2003) models, here- 
after BC03, and other models available in the literature, and indicate systematically 
lower ages and, on average, 60% lower masses for the stellar populations sampled in 



these galaxies. According to Maraston et al. (2006) the source of this discrepancy is pri- 



marily a consequence of the different treatment of the TP-AGB phase in the evolutionary 
models. 



Marigo et al. (2007) have recently concluded new calculations of the TP-AGB evolu- 



tionary phase for stars of different mass and metallicity. The evolution of the stars is now 
computed accounting for the changes in the chemical composition of the envelopes. As a 
consequence of this prescription, the signature of TP-AGB stars around 1 Gyr, i.e the red 
color of the integrated stellar population, becomes more relevant than in previous compu- 



tations. The full implementation of the new stellar libraries and the Marigo et al. (2007) 



TP-AGB evolutionary tracks in population synthesis models is in preparation by Chariot 
& Bruzual (2007, CB07 hereafter). All the BC03 and CB07 SSP models shown in this 



paper have been computed for the Chabrier (2003) IMF and the |Westera et al. (2002) 



stellar library. Fig. 1 shows the time dependence of the fraction of the galaxy luminosity 
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Figure 1. Top frames: time dependence of the fraction of the galaxy luminosity emitted in the 
K-band according to (a) the BC03 and (b) the CB07 SSP models. Different lines correspond to 
different metallicities of the stellar population, as follows: Zq (black solid line), 0.4 x Zq (green 
long-dashed line), 0.2 x Zq (blue short-dashed line), 2.5 x Zq (red dotted l ine). All models s hown 



in this figure were computed for the Padova 1994 tracks assuming the Chabrier (2003) IMF. 
Figures drawn in color are available in the electronic version of this paper. Bottom frames: 1.3 
Gyr Zq isochrone computed using (c) the BC03 and ( d) the CB07 models. The lines representing 
the MS, SGB, RGB, HB, and AGB are identical in both cases, and correspond to the Padova 
1994 tracks. The symbols represent the O-rich (blue squares), the C-rich (black triangles), and 
the superwind (red circles) phases of the TP-AGB. The dotted lines are drawn to guide the eye. 

emitted in the K-band according to the BC03 and CB07 models. From the top frames 
of Fig. 1 it is apparent that the TP-AGB stars in the CB07 models contribute at least a 
factor of two more light in the i-T-band than in the BC03 models, which use a different 
prescription for the TP-AGB evolution (see BC03 for details). At maximum, the TP- 
AGB contributes close to 70% of the K-light in the CB07 model but only 40% in the 
BC03 model. The peak emission in the BC03 model occurs at around 1 Gyr whereas 




Figure 2. Top frames: Percentage of the luminosity of the 1.3 Gyr SSP emitted by TP-AGB 
stars as a function of wavelength (heavy solid- line) according to the CB07 models for (a) 
Z = 0.4 x Zq, and (b) Z = Z@. The contribution of the O-rich TP-AGB stars (red dot- 
ted- line), the C-rich TP-AGB stars (green short-dashed line), and the TP-AGB stars in the 
superwind phase (blue long-dashed line) is also indicated. Figures drawn in color are available 
in the electronic version of this paper. Bottom frames: spectral energy distribution of the 1.3 
Gyr SSP as a function of wavelength (heavy solid-line) according to the CB07 models for (c) 
Z = 0.4 X Zq, and (d) Z = Z Q . The contribution of all the TP-AGB stars (thin solid-line), the 
O-rich TP-AGB stars (red dotted- line), the C-rich TP-AGB stars (green short-dashed line), and 
the TP-AGB stars in the superwind phase (blue long-dashed line) is also indicated. 

in the CB07 model it stays high and close to constant from 0.1 to 1 Gyr. To explore 
the cause of this difference, I plot in the bottom frames of Fig. 1 the 1.3 Gyr Z = Zq 
isochrone for both sets of models. This age corresponds to the spike seen in the Z = Z Q 



line in the top frames of this figure. Whereas the Marigo et al. (2007) tracks used in 
CB07 account for 9 evolutionary stages in the TP-AGB (three in the O-rich phase, three 
in the C-rich phase, and three in the superwind phase), the BC03 models include only 1 
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evolutionary stage on each of these phases. The TP-AGB stars in the CB07 isochrone are 
close to 1 magnitude brighter in K and reach values of V — K several magnitudes redder 
than their counterparts in BC03. The evolutionary rate is such that the total number of 
TP-AGB stars present in the CB07 1.3 Gyr isochrone is 4 times larger than the number 
of these stars present in the BC03 models. The TP-AGB stars represent 0.016% of the 
total number of stars present in this population at this age in the CB07 model, but only 
0.004% in the BC03 model. The net effect of all these factors is the increased contribution 
shown in Fig. 1. 
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Figure 3. Comparison of (a) the B-V and (b) V-K color evolution of the CB07 models (red 
solid line) with the BC03 (blue short-dashed line) and the M05 (black long-dashed line) models. 
Figures drawn in color are available in the electronic version of this paper, (c) and (d): ratio of 
the stellar mass determined from the CB07 and the BC03 models for a given B and Tf-band 
galaxy luminosit y, respectively. Al l models shown in this figure correspond to a Z = Zq SSP 
computed for the |Chabrier (2003)| IMF. 



Fig. 2 shows the contribution of TP-AGB stars to the spectral energy distribution of the 
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Figure 4. Observed spectral energy distribu tion (symbols with error bars) of 7 HUDF galaxies 
taken from Table 1 of |Maraston et al. (20 06) shown together with the best fitting CB07 Z = Zq , 
Chabrier IMF SSP model (red solid-line) . The galaxy HUDF identification number, redshift, age, 
and mass (in units of 10 11 M ) are given inside each frame and listed in Table 1. For clearness, 
after determining the galaxy mass the spectra were forced to share the same vertical scale. 
Figures drawn in color are available in the electronic version of this paper. 



1.3 Gyr SSP as a function of wavelength according to the CB07 models for Z = 0.4 x Zq, 
and Z = Zq. The relative contribution of the O-rich vs. the C-rich TP-AGB stars 
depends critically on the assumed stellar metallicity. Fig. 3 compares the B — V and 
V - K color evolution of the CB07, BC03 and M05 models. In B - V the CB07 and 
BC03 models are identical at all ages. At early and late ages both sets of models have 
the same V — K color, but at intermediate ages the CB07 models are considerably redder 
than the BC03 models. At late ages, the BC03 and CB07 models match very well the 
observations of nearby early-type galaxies, whereas the M05 models are too blue. The 
stellar mass determined from the i\-band luminosity using the CB07 model can be up to 
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Determined from fitting a Z = Zq, Chabrier IMF SSP to the magnitudes 



(2006) 



in Table 1 of Maraston et al. 
2 Values read from Fig. 4 of Maraston et al. (2006)1 assuming E(B — V) 



The galaxy mass Mu is given in units of lO^M 



©• 



50% lower than the mass determined from the BC03 model at the same age. The BC03 
and CB07 models will predict the same galaxy mass if the -B-band luminosity is used. 

Fig. 4 compares the observed spectral energy distribution of 7 HUDF galaxies with 
the best fitting CB07 Z = Z Q , Chabrier IMF SSP model. The fits look quite reasonable 
and can be improved assuming different star formation histories or dust content. The 
galaxy HUDF identification number, redshift, and derived age and mass (in units of 10 11 
M©) are given inside each frame and listed in Table 1. As expected from the previous 
discussion, the mass and age of the stellar population determined from the level of the 
IR flux observed in these galaxies are considerably lower than the values determined by 



Maraston et al. (2006) from the BC03 (and M05) models. The masses derived from the 
CB07 models are 50 to 80% lower than the BC03 masses. The ages derived from the 
CB07 models are, with exceptions, 40 to 60% of those derived from BC03. 
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Discussion 

Kroupa: How could your results on the galaxies (their ages and masses) change if the 
IGIMF is steeper than the standard IMF?. 

Bruzual: If the number of intermediate mass stars is reduced as a consequence of having 
a different IMF, we expect a lower contribution of TP-AGB stars to the galaxy luminosity. 
I would guess that the estimates of the galaxy mass will be somewhat higher than for 
the standard IMF. The age estimates could also show a trend towards larger values. The 
opposite should occur for an increased number of intermediate mass stars in the IMF. 
This is something that can be easily tested with the population synthesis models. 

Cervino: When you use the TP-AGB tracks you compute a mean value, but there will be 
some scatter (TP-AGB stars show variations of a few magnitudes). Have you estimated 
how large this scatter is? How do you manage to integrate isochrones when they are a 
"band" and not a single well defined path? How does the TP-AGB affect the surface 
brightness fluctuations? 



Bruzual: The evolutionary tracks for TP-AGB stars computed by Marigo et al. (2007) 
provide well defined values of (age, log T e ff, Mbol) at all points along the tracks. 
We use their results as published. Presumably the aspects you mention have already 
been taken into consideration by the authors of this paper in their careful comparison 
with observations of real TP-AGB stars. Our isochrones are thus a narrow line and 
not a "band" . At the moment I am working with Rosa Gonzalez-Lopez-Lira in a paper 
where we examine the influence of the new TP-AGB evolutionary tracks on the surface 
brightness fluctuations. 



Maraston: This is not true. There is no scatter. 



Bruzual: Good to hear that. 



Maraston: What differences do you find between the IUE and the HNGSL spectral 
libraries? 



Bruzual: The signal-to-noise-ratio is considerably higher in the HNGSL data than in the 
IUE spectra. As a consequence, the spectral features, absorption lines and discontinuities 
are much sharper and better defined. No doubt, the HNGSL is more appropriate than 
the IUE data to study line strength indices in the UV 

GuSTAFSSON: You mentioned that your high-resolution indices are quite sensitive to the 
degree of degrading due to the stellar velocity dispersion. I suppose that in a real galaxy 
different types of stars may have different velocity dispersion. Is this a worry? 



Bruzual: This is certainly a worry when modeling indices like /200 and ^275 defined 
by Vazdekis & Arimoto (1999, ApJ, 525, 144), which are very sensitive to the stellar 
velocity dispersion. Models shown in Bruzual (2006) assume that all the stars share the 
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same value of the velocity dispersion whereas in reality on should consider the dynamics 
of different stellar groups and use the corresponding a for each kind of stars. 

Han: In your second conclusion, you said that there is room for improvement in UV. 
We know that UV in old stellar populations is from hot subdwarfs (EHB). We have a 
binary model for the formation of hot subdwarfs and the model is well established for 
galactic hot subdwarfs. We applied the model to evolutionary population synthesis study 
and proposed a binary model for the UV-upturn of elliptical galaxies (Poster V.10) 

Bruzual: This is certainly an interesting approach. Binary star evolution is missing from 
most, if not all, stellar population models. The source of the UV in old populations is not 
well established yet. Besides EHB stars in some galaxies there is residual star formation. 
In my talk I referred to the improvement that is needed in the UV spectral libraries, but 
for sure we also need better evolutionary schemes. 



